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Abstract 

A systematic first-principle study is performed to calculate the lattice parameters, electronic 



i-rt '. structure, and thermodynamic properties of UN using the local-density approximation (LDA)-l-C/ 

Q ' and the generalized gradient approximation (GGA)+ U formalisms. To properly describe the strong 

cj : 

correlation in the U 5/ electrons, we optimized the U parameter in calculating the total energy, 
lattice parameters, and bulk modulus at the nonmagnetic (NM), ferromagnetic (FM), and antifer- 
romagnetic (AFM) configurations. Our results show that by choosing the Hubbard U around 2 
l/^ ' eV within the GGA-I- U approach, it is promising to correctly and consistently describe the above 

^—N . mentioned properties of UN. The localization behavior of 5/ electrons is found to be stronger than 



that of UC and our electronic analysis indicates that the effective charge of UN can be represented 
as U^'^^^N^'^^^. As for the thermodynamic study, the phonon dispersion illustrates the stability 
of UN and we further predict the lattice vibration energy, thermal expansion, and specific heat 
by utilizing the quasiharmonic approximation. Our calculated specific heat is well consistent with 
experiments. 
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I. INTRODUCTION 

Uranium nitrides have been extensively studied in experiments in connection with their 
potential applications in the Generation-IV reactors [l]. These reactors raise a number 
of concerns surrounding the issue of nuclear energy. The fission reactions depend on fast 
neutrons, requiring a small core with a high power density and very efficient heat transfer. 
The oxide based fuels are therefore being involved in the ongoing research and development, 
however, the nitride fuels also participate in the competition to become the alternative 
materials for their superior thermal physical properties, such as high melting point, high 
thermal conductivity, and high metal density |2[, as well as the good compatibility with the 
coolant (Na). 

On account of these obvious importances, several studies, such as electronic structure 
optimization |3|, Ij], magnetic properties iSj], point defects ^], and elastic constants 7[, have 
already been conducted for uranium nitride. However, conventional density functional theory 
(DFT) that apply the LDA or GGA underestimates the strong on-site Coulomb repulsion 
of the 5/-electron and, consequently, describes UN as incorrect FM conductor instead of 



a 



the experimentally observed AFM type-I structure [8| at the Neel temperature TAr=53 K. 
Similar problems have been confirmed in studying other electronically correlated materials 
within the pure LDA/GGA scheme. In the present work, we use the LDA/GGA+ U method 



I9 



developed by Dudarev et al. [9| to effectively remedy the failures raised by LDA/GGA in 
describing the strong intra-atomic Coulomb interaction. This method has been successfully 



used to study the correlated problems 10l-ll2| and the total LDA/GGA energy functional is 
of the form 

£^LDA(GGA)+U = ^LDA(GGA) + ^— J^iTrp'^ - Tr(p>'^)], 

where p'^ is the density matrix of / states with spin o", while U and J are the spherically 
averaged screened Coulomb energy and the exchange energy, respectively. 

In this paper, we have systematically calculated the lattice parameters, electronic 
structure, as well as the thermodynamic properties of UN using the above mentioned 
LDA/GGA+ U scheme. We have carefully discussed how these properties are affected by 
the choice of U as well as the choice of exchange-correlation potential. After testing the 
validity of the ground state by choosing U around 2 eV within the GGA+ U approach, we 



performed a series of calculations on the electronic structures, bonding properties, and the 
phonon dispersion. The lattice vibration energy, thermal expansion, and specific heat were 
obtained by utilizing the quasiharmonic approximation (QHA) based on the first-principles 
phonon density of state (DOS). The rest of the paper is organized as follows. The computa- 
tional details of first-principles are briefly introduced in Sec. II. The calculation results are 
presented and discussed in Sec. III. Finally, we give a summary of this work in Sec. IV. 

II. COMPUTATIONAL METHODS 

The DFT total energy calculations were carried out using the Vienna ab initio simulations 



package (VASP) [l3|, [ij] with the projected-augmented-wave (PAW) pseudopotentials [15 1 



and plane waves. The exchange and correlation effects were described within LDA and GGA 
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171]. The uranium 6s^6p^6d'^5f'^7s^ and nitrogen 2s'^2p^ electrons were treated as valence 
electrons. The electron wave function was expanded in plane waves up to a cutoff energy 
of 500 eV and all atoms were fully relaxed until the Hellmann-Feynman forces become less 



than 0.001 eV/A. The Monkhorst-Pack [18| 9x9x9 mesh (75 irreducible k points) was used 
in Brillouin zone (BZ) integration and the corresponding electronic DOS was obtained with 
15x15x15 (120 irreducible k points) A;-point mesh. The strong on-site Coulomb repulsion 
among the localized U-5/ electrons was described by using the formalism developed by 
Dudarev et al. [9|]. In this paper, the Coulomb parameter U was treated as one variable, 
while the parameter J was set to 0.51 eV. Since only the difference between U and J is 
meaningful in Dudarev's approach, therefore, we label them as one single parameter U for 
simplicity. 

III. RESULTS AND DISCUSSIONS 

A. Atomic and electronic structures of UN 

We study UN in its ground-state NaCl-type {Fm3m) structure. In the present 
LDA/GGA-I- U approach, we have considered the NM, FM, and AFM phases for each choice 
of the value of U, and determined the ground-state phase by a subsequent total-energy com- 
parison of these three phases. Compared with the FM and AFM phases, the NM phase is 
not energetically favorable both in the LDA+ U and GGA+ U formalisms. Therefore, the 
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FIG. 1: (Color online) Dependence of the total energies (per formula unit) on U for FM and AFM 

UN. 

results of NM are not presented in the following. The dependence of the total energies (per 
formula unit) for UN in both AFM and FM configurations on U are shown in FigHJ At 
U=0 eV, the ground state of UN is determined to be a FM metal, which is in contrast 
to experiment results. By increasing the amplitude of U, our LDA/GGA+ U calculation 
correctly predicts an AFM metal ground state and the turning value of U is ~1.5 eV and 
~2 eV in GGA+ U and LDA+ U approaches, respectively. In the discussion that follows, 
we, therefore, confine our report to the AFM phase of UN. 

In this paper, the theoretical equilibrium volume Vq, bulk modulus B are obtained by 
fitting the third-order Birch- Murnaghan equation of state (EOS) [19]. Our calculated lattice 
parameter % and B for the cubic unit cells of UN are shown in Fig. [2l For the pure 
DFT calculations (f/=0 eV), both the LDA and GGA methods underestimate the lattice 
parameter with respect to the experimental value. This trend is more evident for LDA 
approach due to its over-binding character. After turning on the Hubbard parameter U, 
the value of ao gradually improves for both LDA and GGA approaches. At around [/=1~2 
eV, the GGA+ U gives ao=4.896~4.926 A, which consists well with the experimental data 



20| of ao=4.888 A. Within LDA+t/, the lattice constant can be satisfied by turning on the 



Hubbard U parameter at around 3 eV. The dependence of bulk modulus i? on t^ is presented 
in Fig. Mjo). It is clear that the LDA+ U results are always higher than that from GGA+ U. 
This is due to above mentioned overbinding effect of the LDA approach. With increasing 
the amplitude of U, the value of B shows a clear declining trend for both schemes. For AFM 
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FIG. 2: (Color online) Dependence of the equilibrium lattice parameter % (a) and the bulk modulus 
B (b) on U. 



phase, the GGA results show that the variety of B is small in the range of f/=2~5 eV . At 
U= 2 eV, the value of B equals to 194.5 GPa, which coincides well with the experimental 
data (194 GPa in Ref. [aj, 200 GPa in Ref. [71, and 206 GPa in Ref. [2^). The LDA 
results always hold higher B values than experimental data till the amplitude of U over 4 
eV. 

On the whole, considering the magnetic configurations, the GGA+ U can give a satis- 
factory prediction of ground-state atomic structures and bulk modulus B by tuning U to 
be near 2 eV for UN. Besides the above effect of LDA/GGA+[/ on the atomic structure 
parameters, in the following discussion we further systematically investigate the electronic 
structures within the two theoretical treatments. The total electronic DOS together with 
the orbital-resolved site-projected DOS (PDOS) of UN are displayed in Fig. [31 Evidently, a 
large degree of U / orbitals can be observed in the valence band near the Fermi level, and 
the conduction band is also strongly marked by / orbitals. Under the SIC-LSD calculations 
in Ref. 23|, the localized and delocalized /-electron configurations were discussed, which 
indicate that the /^ is the energetically favorable configuration. However, it still remains un- 
clear whether a localized, delocalized, or dual localized/delocalized picture can best account 
for the experimentally observed properties of UN. It is only certain that the localization 



in UN is stronger than that in UC 
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231 . For UC, one part of the 5/ electrons transfer 



into the interstitial zone, the other part are expected to be confined to the j=5/2 multiplet, 
and the itineracy of 5/ electrons are evident. Similar to UC, for UN, as the increase of the 
Hubbard parameter [/, the split of j=5/2 and j=7/2 multiplets can also be observed but 
not so clear, and the itineracy of the 5/ electrons still exists in UN. At a typical value of 
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FIG. 3: (Color online) The total DOS for the UN AFM phase computed in the GGA, GGA+U 
(U=2), LDA, and LDA+U (U=4) formalisms. The projected DOSs for the U 5//6d and N 2p 
orbit als are also shown. The Fermi energy level is set at zero. 




FIG. 4: (Color online) Valence charge density of UN (001) plane within GGA+U approach at 
U=2 eV. 

U=2 eV, the conduction band /-electron occupancy is ~2.44 electrons, compared well with 
the 2.2±0.5 electrons measured by Norton et al. J2J] using the photoelectron-spectroscopic 
method. Due to the strong overlap of the U 5/-orbitals near the Fermi energy, the UN phase 
exhibits a clear metallic behavior. 

In order to further analyze the chemical bonding nature of UN, we calculate the effec- 
tive Bader charges [25| in the GGA+[/ formalism with U=2 eV. We adopt 336x336x336 
charge density grids and the spacing between adjacent grid points is 0.011 A. The calcu- 



n 

TABLE I: Bader effective atomic charges of UN. The calculated results using LCAO (Ref. |^]) and 



PW91 (Ref. [26]) methods are also listed for comparison. 



Methods Bader charge UN 

GGA+f7 Qu +1.71 

Qn -1.71 

LCAO Qu +1-58 

Qn -1.58 

PW91 Qu +1-66 

g„ -1.66 

ated valance charges within GGA+ U are listed in Table |T] together with the LCAO (Ref. 

3|) and PW91 (Ref. [26|) results for comparison. Our present Bader analysis gives the 
valency of U^-'^^+N^-^^-, in qualitative agreement with the LCAO (ui-^s+N^-^^-) and PW91 
^1ji.66+]\j-i.66-^ results. This agreement shows that the ionicity in the U-N bond is intrinsic 
and therefore insensitive to the different choices of computation methods. The dominant 
ionic contribution to the U-N bond can also be seen from the total charge density, which is 
plotted in Fig. HI 

B. Phonon dispersion curve of UN 

Through the above discuss on atomic and electronic structures on U, we choose the 
GGA+ U approach with the Hubbard U= 2 eV to calculate the phonon dispersions for UN. 
In calculating the phonon dispersion curves and the phonon density of states, the Hellmann- 
Feynman theorem and the direct method [27| are employed. For the BZ integration, the 
3x3x3 Monkhorst-Pack A;-point mesh is used for the 2x2x2 UN supercell containing 64 
atoms. In order to calculate the Hellmann-Feynman forces, we displace two atoms (one U 
and one N atoms) from their equilibrium positions and the amplitude of all the displacements 
is 0.03 A. Besides, we have calculated the Born effective charges of UN for their critical 
importance in correcting the LO-TO splitting. Because of its high symmetry for UN, the 
off-diagonal elements of the Born effective charge tensor are all zero and the three diagonal 
elements Z^^, Zyy and Zzz are the same. Therefore, only Z^^ is shown here. Our calculated 
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FIG. 5: (Color online) Calculated phonon dispersion curves (left panel) and corresponding DOS 
(right panel) within GGA+ U approach at U=2 eV for AFM UN. 



results for UN are ^^=+1.95 and Z^= — 1.95. The calculated phonon dispersion curves 
along T — X — K — T — L directions is displayed in Fig. O The experimental data from 
Ref. 28| (at T=4.2 K) are also presented for comparison. For NaCl type UN, there are only 
two atoms in its formula unit, therefore, six phonon modes exist in the dispersion relations. 
As shown in Fig. [5l our calculated LA/TA branch is in good agreement with experiment. 
The remarkable splitting between LO and TO at F point can be attributed to the inclusion 
of the Born effective charges in our phonon dispersion calculation. The TO frequency at 
F point is 12.19 THz. This result is well consistent with the available experimental value 
of 12.3 THz at 4.2 K. In addition, the phonon DOS splits into two parts with one part in 
range of 0-4.7 THz where the vibrations of uranium atoms are dominant and another part 
in the domain of 10.5-15 THz where the vibrations mainly come from nitride atoms. This 
evident gap between the optic modes and the acoustic branches is because of the fact that 
the uranium atom is heavier than nitride atom. In the following discussions, the reliability 
of the phonon dispersion calculation will give an accuracy evaluation of the thermodynamic 
properties. 



C. Thermodynamic properties 

To calculate thermodynamical quantities such as the lattice vibration energy, thermal 
expansion, and specific heat, the Helmholtz free energy F in QHA is investigated as follows: 

FiV, T) = EiV) + F,h{y, T) + F,i,{y, T), (1) 

where E{V) stands for the ground state energy, Fph{V,T) is the phonon free energy at a 
given unit cell volume V, and Feie is electron excitation energy. Under QHA, the Fph{V,T) 
can be calculated from phonon DOS by 

Fp^{V,T) = k^T j g(oj)\n[2sinh{^^)]du, (2) 

where u = uj{V) denotes the volume-dependent phonon frequencies, g{u) is the phonon 
density of states, h is the Planck constant, and ks is the Boltzmann constant. Equation ([2]) 
contains some effect of unharmonics since the phonon frequencies have to be derived each 
time at the current crystal volume V. In addition, the specific heat for constant volume Cy 
can be obtained directly as 

Then the specific heat at constant pressure Cp is given by 

Cp-Cv = al{T)B{T)V{T)T, (4) 

where the constant volume thermal expansion ay is defined by ay = ■^(|^)p. The electronic 
excitation effect on the specific heat is accounted by the free-electron Fermi gas model, 
Ce=jTe, where 7 is the electronic specific heat coefficient. For noninteracting electrons, the 
value of 7 is reasonable at low electron temperature (re<3000 K). It is proportional to the 
total density of states N(Ep) at the Fermi level and is given by ^k'^N{Ep). 

In our calculations, the unit cells are expanded and compressed to a set of constant 
volumes in calculations of free energy. Then the equilibrium volume at temperature T is 
obtained by minimizing the free energy. The calculated free energy versus volume curves for a 
number of selected temperatures is plotted in Fig. |6l from which the volume expansion upon 
the temperature increase can be derived. The calculated heat capacity Cp of UN is displayed 



Cv = i^)v = kB I dug{io){^) \ ^7:'-J\,^ . (3) 



in Fig. [3 For comparison, the experimental data from Refs. 29|-|34| and the theoretical 
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FIG. 6: (Color online) Dependence of the Helmholtz free energy F{T,V) on crystal volume at 
various temperatures. The locus of the minimum of the free energy for UN is also presented. 



100 



80 



60 



40 



o 



20 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



This work (lattice) 

This worl< (lattice + electronic) 




Counsell etal. (1966) 
Affortit(1969) 
Getting etal. (1972) 
Tagawa (1974) 
Matsuietal. (1987) 
Chevalier etal. (2000) 
I I I I I I I I I I I I I I 



+ Week etal. 

(2007 lattice) 
Week et al 

(lattice+electronic)_ 



500 1000 1500 2000 2500 3000 

T(K) 



FIG. 7: (Color online) Heat capacity of UN calculated within GGA+ U at U=2 eV with and with- 



out considering the contribution of electrons. Experimental data from Refs 



results from Refs. 



35|, 



29l434| and theoretical 



36( 1 are also displayed for comparison. 



results by Chevalier et al. (N-U modelling calculation) [35[ and Week et al. [36| (all-electron 
calculation) are also presented. As shown in Fig. U\ the calculated thermodynamic functions 
with no electronic excitation contribution at low temperature are in good agreement with the 
experimental data up to around 200 K. However, in the high temperature domain (r> 200 
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K), our calculated results with only lattice vibration included evidently underestimate the Cp 
compared with experimental values. This kind of underestimation in the high temperature 
domain has also been observed in the all-electron calculation 36|]. Therefore, one needs to 
take into account the conduction electrons contributions to the Cp for metallic material UN. 
Our estimated value for the electronic specific heat coefficient 7 is equal to 26.7 mJ K^^ 
mol~^. Although this value is somewhat lower than the experimental value of 49.6 mJ K~^ 
mol~^ |37|, it is apparently larger than the value of 3.1 mJ K ~^ mol~^ for its isostructural 
analog ThN [38|. This can be contributed to the unfilled 5/ electrons for UN. The specific 
heat capacity including electronic contribution is also displayed in Fig. [71 One can see that 
the Cp with electronic corrections is largely enhanced, in agreement with the experimental 
values within a broad temperature domain. Therefore, the partially itinerant U 5/ electrons 
not only play a great role in electronic properties, but also have considerable infiuence on 
the thermodynamic performances of the intermetallic compounds such as the present UN. 

IV. CONCLUSION 

In summary, we perform systematic first-principles calculations on the structural, elec- 
tronic, and thermal properties of UN using the LDA/GGA-I- U method. With the Hubbard 
U correction, the antiferromagnetic nature of UN is successfully predicted. The atomic 
structure, including lattice parameters and bulk modulus can be reasonably given, com- 
pared with corresponding experimental values. The calculated electronic density of states 
shows the important role that the 5/ electrons play in the conduction band as well as in 
the valence band. By choosing U=2 eV within GGA, the phonon dispersions and phonon 
density of states can be reasonably derived with regard to the experimental data. Using the 
reliable phonon spectrum, the specific heat for constant pressure C p of UN including both 
lattice and conduction electron contributions are calculated, the results of which are in good 
agreement with experimental data. We expect that our calculated results will be useful for 
the application of uranium nitride in the Generation-IV reactor and nuclear industry. 
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